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Abstract 

I review recent results and unsolved problems concerning the hard-core 
lattice gas and the g-coloring model (antiferromagnetic Potts model at zero 
temperature). For each model, I consider its equilibrium properties (unique- 
ness/nonuniqueness of the infinite- volume Gibbs measure, complex zeros of the 
partition function) and the dynamics of local and nonlocal Monte Carlo algo- 
rithms (ergodicity, rapid mixing, mixing at complex fugacity). These problems 
touch on mathematical physics, probability, combinatorics and theoretical com- 
puter science. 

[Pour M. Toubon] Je passe en revue des resultats recents et des problemes 
non resolus concernant le gaz sur reseau avec exclusion et le modele de q- 
coloriage (modele de Potts antiferromagnetique a temperature nulle). Pour 
chacun des deux modeles, je considere ses proprietes d'equilibre (unicite/non- 
unicite de la mesure de Gibbs en volume infini, zeros complexes de la fonction de 
partition) et la dynamique des algorithmes Monte Carlo locaux et non-locaux 
(ergodicite, melange rapide, melange a fugacite complexe). Ces problemes 
touchent a la physique mathematique, a la probabilite, a la combinatoire et 
a I'informatique theorique. 
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1 Introduction 



In this paper I propose to review some recent results and to list some open prob- 
lems concerning a pair of statistical-mechanical lattice models (defined precisely in 
Sections and H): 

• Hard-core lattice gas with nearest-neighbor exclusion and fugacity w 

[= "independent- set model" in graph-theory language] 

• g-coloring of the vertices of a graph 

(These are zero-temperature models. More generally, one can study the hard-core 
lattice gas with a soft nearest-neighbor repulsion, and the antiferromagnetic Potts 
model at nonzero temperature.) 

What these two models have in common is that they are "antiferromagnetic": 
nearest neighbors prefer to be (in the g-coloring case have to be) in different states. 
Unlike ferromagnetic models, for which (roughly speaking) only one type of ordered 
phase is possible, in antiferromagnetic models the possible types of order can depend 
delicately on the parameters (e.g. w or q) and on the nature of the lattice (e.g. bipartite 
or not). This motivates studying these models on an arbitrary graph G = {V,E) — 
and not just on a regular lattice — in order to investigate which features of G correlate 
with which features of the phase transition. Here G will be either finite or countably 
infinite, depending on the question being investigated. 

There is one other motivation for studying these models on arbitrary graphs: The 
hard-core lattice gas on a general graph G (with different fugacities G C assigned 
to different vertices) is in fact the universal statistical-mechanical model in the sense 
that any statistical-mechanical model living on a vertex set Vq can be mapped onto a 
gas of nonoverlapping "polymers" on Vq, i.e. a hard-core lattice gas on the intersection 
graph of Vq Section 5.7]. 

For each of these models, I will discuss the following questions: 

• Equilibrium statistical mechanics 

— Uniqueness/nonuniqueness/properties of Gibbs measures in infinite vol- 
ume 

— Complex zeros (in w or q) of the partition function in finite volume 

• Monte Carlo algorithms (single-site, local and nonlocal) 

— Ergodicity (this is sometimes nontrivial at zero temperature) 

— Rapid mixing 

— Mixing at complex i;; or g (a new and perhaps crazy idea) 

I want to stress that I'm not an expert on any of these things, so some of my questions 
may be naive and some of the "open problems" that I pose may already be solved. 
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2 Lattice Gas: Equilibrium Properties 



Let G = (V, E) be a finite graph, and let w = {wx}x£V be a set of (possibly com- 
plex) fugacities. The partition function of the hard-core lattice gas (= independent- set 
polynomial) on G is 

ACV xgA 
A independent 

(Recall that A (1 V is called an independent set if it does not contain any pair of 
adjacent vertices.) When the fugacities are nonnegative, an equilibrium probability 
distribution on the state space 2^ can be defined by 

/i(A) = Z(w;)"^JJw^ X x(^ independent) . (2.2) 

xGA 

When G is countably infinite and the fugacities are nonnegative, Gibbs measures can 
then be defined by the usual Dobrushin-Lanford-Ruelle prescription . 

At small w we expect to have uniqueness of the Gibbs measure, exponential decay 
of correlations, analyticity, etc. Indeed, the Dobrushin uniqueness theorem [^, ^ 
easily implies the uniqueness of the infinite-volume Gibbs measure, and the expo- 
nential decay of correlations in this unique Gibbs measure, whenever < < 
{1 — e) / {dx — I + e) for all vertices x [here e > 0, and dx denotes the number of vertices 
adjacent to x]. In particular, when all vertices are assigned the same fugacity w, this 
happens whenever 0<w<l/(A — 1) [here A = maxx^v dx is the maximum degree 
of G]. 

A slightly better bound |Q| can be obtained using the disagreement-percolation 
method: if the site percolation model on G with occupation probabilities px = 
Wx/{1 + Wx) has (with probability 1) no infinite occupied cluster, then the hard-core 
lattice gas with fugacities w has a unique Gibbs measure; moreover, the two-point 
covariances in the lattice gas are bounded above by the corresponding connection 
probabilities in the percolation model. In particular, if < < (1 — e)/{dx — 1) for 
all vertices x, then standard arguments (comparison to a branching process) show that 
the percolation model has (w.p. 1) no infinite occupied cluster and has exponential 
decay of connectivities. So we prove uniqueness and exponential decay for the lattice 
gas whenever < < (1 — e)/ {dx — 2 + e) for all vertices x, or < w < 1/(A — 2) 
in case the fugacities are all equal. 

It is natural to ask whether this result is sharp. The answer is almost certainly 
no: Vigoda [^, ^ improves 1/(A — 2) to 2/(A — 2), at least for regular lattices, as a 



byproduct of proving rapid mixing of Glauber dynamics (see Section ^ below) .[| But 
I suspect that this result is not sharp either, and that the optimal result is: 



^ For some lattices, the disagreement-percolation bound is better than 2/(A — 2). For example, 
the triangular-lattice site-percolation model has Pc = 1/2 Js^ , which implies uniqueness of the 
triangular-lattice-gas Gibbs measure for for w < 1; but 2/(A — 2) = 1/2. 
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Conjecture 2.1 For any countably infinite graph G of maximum degree A, the hard- 
core lattice gas on G has a unique Gibbs measure whenever < < {1 — e)x 
(A — l)^~-^/(A — 2)^ [~ e/A for large A] for all vertices x, for some e > 0. [Perhaps 
this is true even with e = 0./ 

Conjecture is motivated by the fact that if = (A — 1)^^^/(A — 2)^ is the critical 
point for the complete rooted tree with branching factor r = A — 1 ||5^, |^ ^ |65| . 
It seems to be an open problem to prove Conjecture |2.1| even when G is a non-regular 
tree and the are all equal. 

On a finite graph G, the free energy log Z{w) and the correlation functions (which 
are derivatives of log Z) are analytic functions of w on any simply connected domain 
D G C where Z{w) is nonvanishing. Moreover, by a standard Vitali argument |^ 



pp. 417-418], this analyticity can be carried over to the infinite- volume limit whenever 
this limit exists for real w (e.g. for amenable transitive graphs). So it is of great 
interest to find sufficient conditions for the nonvanishing of Z{w) [which is trivial 
when w > 0, but not when w is negative or complex]. Shortly before his death, 
Dobrushin |T^, [1^ proved a beautiful theorem on the absence of zeros of Z{w) at 
small complex fugacity: 

Theorem 2.2 (Dobrushin [|T8|, |19|| ) Let{cx}xGV be an arbitrary set of nonnegative 
numbers, and define 



Rx = (1-e "^) expl-^cj 



(2.3) 



where ?/ ~ x denotes that y is adjacent to x. Then Z{w) is nonvanishing in the closed 
poly disc \wx\ < Rx- 

The proof is an extraordinarily simple 3-line induction on the number of vertices in 
G: see |l^, 19|, or see [^, Section 3] for an equally simple proof of a slightly stronger 
result. 

Corollary 2.3 IfG has maximum degree A, then Z{w) is nonvanishing in the closed 
polydisc \wx\ < A^/(A + 1)^+^ [~ l/(eA) for large A]. 



Remarks. 1. Kotecky and Preiss proved a result slightly weaker than Theo- 



rem |2]^, in which 1 — e '^'^ is replaced by CxC 

2. Corollary |2.3| is close to best possible: for the complete rooted tree with 



branching factor r = A — 1 and depth n, Z{w) has negative real zeros that tend to 
w = —{A — 1)^~^/A^ as n — > oo Indeed, by an inductive argument closely 

related to (but more subtle than) Dobrushin's, Shearer |Q improves Corollary p.3| 



to the optimal radius (A — 1)^~^/A^. Alex Scott and I are investigating analogous 
improvements of Theorem |2.2| . Note that this optimal radius is a factor ^ e smaller 
than the interval at real positive w that is obtained by the Dobrushin uniqueness 
theorem, and a factor ^ smaller than the interval proposed in Conjecture p.l| . 
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3. For any lattice gas with repulsive interactions — and in particular for one with 
hard-core exclusions — the Mayer expansion log Z{w) = '^c^w"^ has alternating 
signs: (— l)'"'~^Cn > 0. It follows that the zero of Z{w) closest to the origin always 
lies on the negative real axis. 



4. Theorem p.2| (which can be generalized to soft-core repulsive interactions pi 
provides an extraordinarily simple proof of the convergence of the Mayer expansion 
for such lattice gases. Recall that the usual approach to proving convergence of the 



Mayer expansion [46, 57, 16, 11, 59, 12, p4] is to exphcitly bound the expansion coeffi- 



cients; this requires some rather nontrivial combinatorics (for example, an inequality 



of Penrose together with the counting of trees). Once this is done, an immediate 
consequence is that Z is nonvanishing in any polydisc where the series for logZ is 
convergent. Dobrushin's brilliant idea was to prove these two results in the opposite 
order: first one proves, by an elementary induction on the cardinality of the state 
space, that Z is nonvanishing in some specified polydisc; it then follows immediately 
that log Z is analytic in that polydisc, and hence that its Taylor series is convergent 
there. It is an interesting open question to know whether this approach can be made 
to work without the assumption of hard-core self-repulsion. 

What if we ask about regions of if-space other than polydiscs centered at w = 0? 

Question 2.4 What is the largest complex domain containing w = on which Z{w) 
is zero-free for all graphs of maximum degree A? [This question has two versions, 
depending on whether the are assumed equal or unequal. In the latter version, one 
seeks domains D for which Z{w) ^ whenever Wx & D for all x; hut there need not 
exist a unique maximal such domain.] 

I conjecture that there is a complex domain Z^a containing at least the interval 
< w < 1/(A — 1) of the real axis — and possibly even the interval Q < w < 
(A — 1)^^^/ (A — 2)^ — on which Z{w) is zero-free for all graphs of maximum degree 
A. See also Section |^ regarding mixing at complex fugacity. 

Another question is: Can we get sharper constraints on the location of the complex 
zeros of Z{w) for restricted classes of graphs Gl A classic result of this kind is due 
to Heilmann and Lieb (see also p7| , |5I|): if G is a line graph., then all the zeros 
of Z{w) are negative real. [The lattice gas on a line graph L{H) is identical to the 
matching polynomial (monomer-dimer-model partition function) on H.] Now, it is 
well known that G is a line graph if and only if nine particular graphs Gi, . . . , Gg do 
not appear as induced subgraphs of G [^, pp. 73-77]. It turns out that for eight of 
these "forbidden induced subgraphs", all the zeros of Z{w) are negative real; the only 

one with non-real zeros is the claw Ki^^ = * This suggests the conjecture: 



Conjecture 2.5 (Hamidoune ||28| , Stanley [|66|| ) // G is claw-free (i.e. has no 
induced subgraph Ki-^), then all the zeros of Z{w) are negative real. 



5 



As vague evidence in favor of this conjecture, let us mention the following: (a) It 
implies the log concavity of the counts of independent sets, which is indeed true 
for claw-free graphs [Q. (b) It is implied by a plausible conjecture concerning a 
special symmetric function Xq associated to G |^ Corollary 2.10]. (c) The recent 
disagreement-percolation proof of complete analyticity for the monomer-dimer model 
generalizes immediately to the hard-core lattice gas on any claw-free graph 
For possibly useful background on claw-free graphs, see see also ||2^. For 
interesting results along vaguely related lines, see |51|, IST 



Finally, one may ask: What happens when the fugacity w is not small? In par- 
ticular, one would like to understand the uniqueness or nonuniqueness of the infinite- 
volume Gibbs measure as a function of w, the nature of the large-w Gibbs measures, 
and whether the phase transition(s) is/are first-order or second-order. On a bipartite 
graph with symmetry between the two sublattices, the mechanism driving the phase 
transition is clear: at large w, one sublattice becomes preferentially occupied and the 
other preferentially vacant. But this does not necessarily exhaust the possible Gibbs 
measures: for example, on the 3-regular tree there also exist Gibbs measures in which 
the density is not constant on each sublattice 0. Moreover, the nonuniqueness of 
the Gibbs measure is not necessarily monotone in i;; |p. For general (non-bipartite) 
graphs, the situation is even worse, and the nature of the possible ordered phases 
seems poorly understood at present. 



3 g-Coloring Model: Equilibrium Properties 

Let G = {V, E) be a finite graph and let g be a positive integer. A map a: V 
{1,2, . . . ,q} is called a proper q-coloring of G in case a{x) ^ for all pairs of 
adjacent vertices x ~ We define /i to be the uniform probability distribution 
on the set of proper g-colorings of G (whenever this set is nonempty). When G 
is countably infinite, Gibbs measures can then be defined by the usual Dobrushin- 



Lanford-Ruelle prescription |26[. The g-coloring problem is the zero-temperature 
limit of the antiferromagnetic g-state Potts model |7^, [73[| . 

When q is large compared to the maximum degree A, the conditional probability 
distribution of a{x) depends only weakly on {o"(?/)}j;^x, so we expect to have unique- 
ness of the infinite- volume Gibbs measure and exponential decay of correlations. This 
is in fact the case: 



Theorem 3.1 (Kotecky [g^], Salas and Sokal []5|]) If q > 2A, the hypotheses 
of the Dohrushin uniqueness theorem hold. Consequently, there is a unique infinite- 
volume Gibbs measure, and it has exponential decay of correlations. 

More generally, for antiferromagnetic Potts models with q > 2 A, there is a unique 
infinite-volume Gibbs measure, with exponential decay of correlations, uniformly 
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down to and including zero temperature pSl. Otherwise put, zero temperature be- 



longs to the high-temperature regime! The reason, of course, is the large ground-state 
entropy. 

Salas and Sokal |^ have some improvements of this result for specific regular lat- 
tices (square, triangular, hexagonal, Kagome), based on decimation and a computer- 
assisted proof. But their results are far from sharp. 



Problem 3.2 Prove uniqueness of the infinite-volume Gibbs measure and exponential 
decay of correlations for q = 3 on the hexagonal lattice, q = 4,5,6 on the square lattice, 
etc. 



For the g = 4 square-lattice model, perhaps one could use the Ashkin- Teller repre- 
sentation plus antipercolation ideas. 

It is natural to ask whether the borderline 2 A is sharp for some (perhaps irregular) 
lattice. The answer is most likely no: Vigoda |^ improves 2 A to yA, at least 



for regular lattices, as a byproduct of proving rapid mixing of a particular local 
dynamics (see Section |] below). But I suspect that this result is not sharp either. 
Even for the A-regular tree, the correct borderline is not known; all that is known 
is that there is a nonunique Gibbs measure whenever q < A ^ and a unique 

Gibbs measure whenever g > a certain qdA) that satisfies qdA) < |A — 1 and 
limsup^^Qo gc(^)/^ < 1.6296 0. For trees the Gibbs measure may well be unique 
as soon as g > A; and it is at least conceivable that this is the sharp borderline for 
general graphs. 



By the Fortuin-Kasteleyn representation |3^, |2J] , the g-state Potts model can be 
analytically continued to complex g; indeed, on any finite graph the Potts-model par- 
tition function is a polynomial in g (see e.g. |]6l| . Section 1] for a review). Specializing 
to the g-coloring problem, we conclude that Pcil) = # (proper g-colorings of G) is 
the restriction to Z_|_ of a polynomial in g, which is called the chromatic polynomial of 



G [|^, [T^. It therefore makes sense to study the real or complex zeros of PgIq)- 



Here are some recent results: 



Theorem 3.3 (Sokal |]6T|]) 

(a) For a graph G of maximum degree A, all the (real or complex) zeros of Pciq) He 
in the disc \q\ < G{A), where G{A) is the solution of a certain transcendental 
equation and satisfies C(A) < 7.963907A. 

(b) For a graph G of second-largest degree A, all the (real or complex) zeros of 
Pg{i) He in the disc \q\ < G{A) + 1. 

The proof is based on writing the Fortuin-Kasteleyn representation of Pg{(i) as a 
gas of nonoverlapping polymers P with fugacity w{P) ~ q~^^^^~^\ When |g| is large 
enough, one can apply Theorem to this polymer gas, and conclude that Z ^ 0. 
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Question 3.4 

(a) What is the sharp hound Copti.^) ? Might 2A work? Is Copt(3) = 3 ? What about 
regions of the q-plane other than discs centered at q = 0? 

(b) Can "second-largest degree" be improved to "max flow"? (See 
for details.) 

(c) Can one obtain a sublinear bound for suitable subclasses of graphs? 

— For theta graphs Qsi,...,sp [consisting of a pair of endvertices connected by 
p internally disjoint paths of lengths Si, . . . , Sp] the answer is yes; there is 
a bound of order p/ log p / |7^/ . 

— For series-parallel graphs, I conjecture that the answer is yes (again p/ logp) 
for maximum degree, and no for second-largest degree or max flow; I am 
currently trying to prove these claims. 

— For planar graphs, it is an open question. 

Let C be the (countably infinite) set consisting of tlie cliromatic roots of all finite 
graphs, and let C be its closure. Jason Brown (unpublished) asked whether C has 
nonzero two-dimensional Lebesgue measure. The answer is yes, in the strongest 
possible sense: C is the whole complex plane! More precisely, one has: 

Theorem 3.5 (Sokal ||63|| ) For theta graphs B*^*'^^ [= p parallel- connected chains 
each consisting of s edges in series], the roots of their chromatic polynomials, taken 
together, are dense in the entire complex q-plane with the possible exception of the 
disc |g — 1| < 1. 

For a brief discussion of the intuition behind the proof, see [p^] . 

Corollary 3.6 (Sokal [|63|| ) There is a countably infinite family of (not-necessarily- 
planar) graphs whose chromatic zeros are, taken together, dense in the entire complex 
q-plane. 

Indeed, it suffices to consider the union of the two families B^^'P^ and O^''"^ + K2 and 
to recall that PG+K„{q) = q{q - I) ■ ■ ■ {q - n + 1) Pciq - n)^ 

^ Here K2 = * • is the complete graph on two vertices. The join of two graphs, denoted 

Gi + G2, is obtained from the disjoint union Gi U G2 by adjoining an edge between each pair of 
vertices a; G Gi and y G G2. 



61, Section 7] 
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4 Lattice Gas: Dynamics of Algorithms 



The single-site heat-bath (Glauber) dynamics for the hard-core lattice gas on a 
finite graph G = (V, E) is defined as follows: Choose uniformly at random a vertex 
X & V; make x occupied with probability Wx/{1 + Wx) if all the neighbors of x are 
vacant, and make x vacant otherwise. This is easily seen to define an irreducible 
Markov chain satisfying detailed balance with respect to the Gibbs distribution /x. 
We want to know: Under what conditions is this Markov chain rapidly mixing, in 
the sense that its mixing timeQ is bounded by a (hopefully low-order) poljTiomial in 
n (the number of vertices in G')?0 

The standard proof of the Dobrushin uniqueness theorem (see e.g. [^, Section 5.1]) 
shows also the O(nlogn) mixing of the Glauber dynamics whenever the Dobrushin 
hypothesis holds. In particular, this holds for the hard-core lattice gas whenever 
< w < 1/(A — 1). This result can alternatively be derived by path-coupling 
arguments, which also show 0(n^ logn) mixing at w = 1/(A — 1) |^ . 



These results are, however, not optimal. By more intricate coupling arguments. 



Luby and Vigoda [^, 37, 68] have recently proven O(nlogn) mixing of the Glauber 
dynamics for < w < 2/(A — 2), and 0{v?) mixing at w = 2/(A — 2).[| 

The mixing properties of the Glauber (or any local) dynamics are closely related 
to the model's equilibrium properties. Indeed, consider a model defined on the regular 
lattice Z'^, and suppose that there exists a local dynamics with mixing time 0(L'^+^~'^) 
on a cube of side L, uniformly in the boundary condition. Then van den Berg 
shows that there is a unique infinite-volume Gibbs measure. Conversely, suitable 
conditions of "weak dependence on boundary conditions" — which are stronger than 
uniqueness of the Gibbs measure, but which are implied by most uniqueness proofs 
— imply O(L'^logL) mixing of the Glauber dynamics, uniformly in the boundary 
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We use here the computer scientists' definition of mixing time: Let 



y 

be the distance to stationarity after time i, starting in state x; and define the (worst-case) mixing 
time 

T(e) = maxmin{i: /^^{t') < e for aU t' >t} . 

X 

Other quantities of interest are the autocorrelation times Ti^t^A and Tcxp defined in ]60[ |; in particular, 
Texp is essentially the inverse spectral gap. [More precisely, exp(— 1/roxp) — spr(P), where spr(P) is 
the spectral radius of P acting on the quotient space modulo constant functions.] 

^ One can also consider more general local dynamics, such as the single-edge heat-bath algorithm 
m, |l| or the Dyer-GreenhiU chain 

^ Dyer and Greenhill j2l] earlier proved O(nlogn) mixing when < w < 2/(A — 2), and 
O(n^logri) mixing when w = 2/(A — 2), for a slightly more complicated dynamics that acts on a 
vertex v and its neighbors. By a comparison argument (see also p^), they then deduced 0{ii? logn) 
mixing when < u; < 2/(A — 2), and 0{n^\ogn) mixing when w — 2/(A — 2), for the Glauber 
dynamics. But the Luby- Vigoda direct analysis appears to give better bounds. 
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condition • 



I'd like to raise a general question: What is the relation between proofs of rapid 
mixing by coupling methods and by the Dobrushin uniqueness method? Otherwise 
put, suppose we try to translate the coupling proof of rapid mixing from probabilistic 
to analytic language: we conclude that the transition matrix P (acting on the quotient 
space modulo constant functions) has spectral radius < /3 < 1, hence has norm < 
(3 + e <1 in some norm or other — but which norm? I can give an answer in the usual 
situation where the coupling is Markovian, with a transition matrix Q{x,y —>■ x',y') 

satisfying J^y'Qi^^V ~^ ^'^u') = for all y and J^x'Qi^^V ~^ ^'^u') = 

P{y — > y') for all x, and where there is a Lyapunov function $(x,?/) > satisfying 

E(<l>{Xt+,,Yt+,)\Xt,Yt^ < mXt,Yt). (4.1) 
Then a simple calculation shows that 

||P/||Lip($) < , (4.2) 

where the Lipschitz seminorm is defined by 

||/||up(.) - -P $(,^,) • (4-3) 

This approach could perhaps be useful in trying to understand the circumstances 
under which the coupling method can or cannot give sharp bounds ||14|| . 

My reasons for wanting to translate a beautiful probabilistic proof into ugly ana- 
lytic language are linked to my next topic, which is the mixing properties of Glauber 
dynamics at complex fugacity w. (I know this sounds crazy, but please bear with 
me.) Since the elements of the Glauber transition matrix P are rational functions of 
w, we can certainly extend them to complex w. The matrix P continues to satisfy 
PI = 1, but its rows P{x, ■ ) are complex measures with total variation norm > 1 
whenever w is negative or complex. Nevertheless, it is easy to prove the following 
result: 

Theorem 4.1 (Sokal ||65|| ) Let G be a finite graph, and suppose that at fugacity 
w e C the Glauber transition matrix P has spectral radius < 1 ( acting on the quotient 
space modulo constant functions). Then Z{w) ^ 0. 

This theorem thus provides a technique for proving the nonvanishing of Z{w) that is 
alternative to Theorem |2.2| , and is inspired in part by old work of Israel |^ . 

The problem is: How to verify that spr(P) < 1? One might hope to do this 
adapting the standard proof of the Dobrushin uniqueness theorem (see e.g. [§3); but 



surprisingly, this proof falls apart immediately as soon as w is even slightly negative 
or complex, if one employs (as usual) the total oscillation seminorm. (The proof fails 
in the strangest way: "dusting" a site i can cause the "dirt" on distant sites j to 
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be amplified exponentially!) But this does not mean that the Dobrushin uniqueness 
method is failing; it may mean only that we are not using the right seminorm. Indeed, 
for any finite graph G we clearly have spr(P) < 1 for all w G [0, oo), so by continuity 
this must hold also in some (G-dependent) complex neighborhood of [0, oo). Which 
one? In particular, can we find a complex domain in which spr(P) < 1 holds for all 
graphs G of maximum degree A? 

Conjecture 4.2 There exists a complex domain Da containing the interval 
[0, 1/(A — 1) ) such that spr(P) < 1 holds for all graphs G of maximum degree A 
whenever w G D^- [As with Question |^.^ , there are two versions of this conjecture, 
depending on whether or not the are assumed equal] 

Indeed, it is not out of the question that 1/(A — 1) can be replaced here by 2/(A — 2) 
or even by (A - 1)^"V(A - 2)^. Unfortunately, what I am thus far able to prove is 
much weaker, namely: 

Theorem 4.3 (Sokal ||65|| ) Let {cx}x£V be an arbitrary set of positive numbers, and 
define 

= (e^--l) expl -J^cJ (4.4) 

where y ^ x denotes that y is adjacent to x. Then spr(P) < 1 holds in the open 
polydisc \wx/{l + Wx)\ < Rx- In particular, Z{w) is nonvanishing there. 

Corollary 4.4 (Sokal [|65|| ) // G has maximum degree A, then spr(P) < 1 and 
Z{w) 7^ hold in the open polydisc 

<^^. (4.5) 



1 + Wx 



Theorem |4.3| is strikingly similar to Theorem |2.2| , but is strictly stronger [since Pj./(1 + 



Rx) > Rx]- Corollary is very close to Shearer's improvement of Corollary p.3| , but 
is slightly weaker (resp. stronger) than Shearer's result when w < (resp. w > 0); 
it is thus almost but not quite sharp for w < 0. The proof of Theorem uses an 
exponentially weighted seminorm on "Fourier coefficients" that is completely different 
from the total oscillation seminorm: Let = 0, 1 be the occupation variable at site 

X ^ V] write /(r/) = ^ ax Yl ~ Vx) and define 

xcv xex 

lll/lllc = E e^p( E^-J l^^l (4.6) 

xcv \ 2;eX / 

Xj^d) 

for each nonnegative vector c = {cx}x€V- The question is: How to interpolate 
smoothly between these two radically different seminorms, which work in different 
regions of w-space, in order to prove spr(P) < 1 for a domain containing both the 
disc (U) and the interval [0, 1/(A - 1) )? 
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5 Q'-Coloring Model: Dynamics of Algorithms 



The single-site heat-bath (Glauber) dynamics for the g-coloring model on a finite 
graph G = {V, E) is defined as follows: Choose uniformly at random a vertex x ^ V, 
and give a new value (independent of the old one) from the uniform distribution 
over all colors different from {o'{y)}y^x- More generally, one can consider the single- 
edge heat-bath dynamics in which a pair of adjacent sites are simultaneously recolored 
conditional on their neighbors pOf, or a yet more complicated local dynamics (e.g. 
1131 , |69[]). Finally, one can consider the Wang-Swendsen-Kotecky (WSK) nonlocal 
cluster dynamics [^0|, |7l|, which is defined as follows: Choose uniformly at random 
a pair of distinct colors G {l,...,g}; let Gais be the induced subgraph of G 
consisting of sites x for which (t{x) = a or /3; then, independently on each connected 
component of G^fB, with probability | either leave that component alone or else 
interchange colors a and /5 on it. 

All these Markov chains are easily seen to satisfy detailed balance with respect to 
the uniform distribution over proper g-colorings. However, since we are at zero tem- 
perature, ergodicity (= irreducibility) is already a nontrivial question. The following 
results are known: 



Single-site dynamics is ergodic for g > A -|- 2 |^3[, and can be nonergodic for 
g = A -|- 1 (consider G = -R'a+i) 



Single-edge dynamics is ergodic for g > A -|- 1 | 
q = A (consider G = -R'a+i minus a single edge) 



], and can be nonergodic for 



WSK dynamics is ergodic for g > A + 1 



WSK dynamics is ergodic for all g if G is bipartite [|T^, ^ 



WSK dynamics is ergodic for g = 4 on subsets of the triangular lattice with free 



boundary conditions |^ 



WSK dynamics is nonergodic for g = 3 on 3m x 3n square lattices with periodic 



boundary conditions, whenever m, n are relatively prime 



WSK dynamics is nonergodic for g = 4 on the 6x6 triangular lattice with 
periodic boundary conditions, and for g = 3 on the 3x3 Kagome lattice with 
periodic boundary conditions^ []54| 

Whenever g < A/2 and g is sufficiently large, there exists a graph of maximum 
degree A and chromatic number 0(g/logg) on which the WSK dynamics for 



g-colorings is nonergodic 43 



^ By a "3 X 3 Kagome lattice" I mean one that is obtained from a 3 x 3 triangular lattice (= 
the Bravais lattice of the Kagome) by placing 3 Kagome sites on each triangular site. It thus has 
3 X 3 X 3 = 27 sites. 
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All other cases are open (to the best of my knowledge). In particular, we would like 
to know for which values of n (if any) the WSK dynamics is ergodic for g = 4 on an 
n X n periodic triangular lattice, or for g = 3 on an n x n Kagome lattice.[] 

When the algorithm is ergodic, the next step is to investigate its mixing rate. 
Here is what is known: 

(a) For the single-site heat-bath dynamics, the Dobrushin uniqueness argument 
proves 0{nlogn) mixing whenever q > 2 A [Q. This result can alternatively be 
derived by coupling arguments, which also show 0{n^) mixing at g = 2 A [^3|, pO| . 
More generally, these results hold for the antiferromagnetic Potts model at inverse 
temperature /? in case q > 2A(1 — e~^) [Q. 

(b) For the single-edge heat-bath dynamics, a path-coupling argument shows 
O {n log n) mixing whenever q > 2 A . 

(c) For a local algorithm involving simultaneous recoloring of a site and all its 
neighbors, a path-coupling argument shows O(nlogn) mixing for q = 5 when A = 3, 
and for q = 7 when G is a 4-regular triangle- free graph ||T^ . This shows that the 2 A 



bound can be beaten. By a comparison argument, one also obtains rapid mixing of the 
single-site heat-bath dynamics, but the bounds are poor [0{n^logn) and 0(n''logn), 
respectively]. 

(d) For a local algorithm involving WSK moves on components of size < 6 only, a 
path-coupling argument shows 0(?2logn) mixing for q > ^A, and polynomial-time 
mixing at g = ^A ||6^. By a comparison argument, one also obtains 0(n^ log n) 
mixing for g > ^A, and polynomial-time mixing at g = yA, for the single-site 
heat-bath dynamics and for the single-cluster variant of the full WSK dynamics 



(e) For the full WSK algorithm, numerical data strongly suggest that there is 
constant-time mixing (i.e. the complete absence of critical slowing-down) for the g = 3 
model on n X n periodic square lattices (at least for n even), uniformly down to and 
including the zero-temperature critical point It would be wonderful to be 

able to prove something! 



(f) Luczak and Vigoda construct a family F„ of planar (hence 4-colorable!) 
graphs on which the WSK algorithm is torpidly mixing [i.e. has mixing time > exp(n'^) 
for some 5 > 0] for arbitrarily large g [indeed, for any g = 0{n^~''^)]. They also show 
that for each pair (g. A) satisfying 3 < g < A/(201og A), there exists a family Gg,A,n 
of bipartite graphs of maximum degree A on which the WSK algorithm is torpidly 
mixing. 
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